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OO . First results from lattice QCD revealing the chiral nonanalytic behavior of nucleon and 
A baryon magnetic moments are presented. Numerical simulations in the light quark 
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O ! mass regime employing the nonperturbatively C(a)-improved conserved vector current 
2 I are enabled via FLIC fermions. Quenched chiral perturbation theory for the nucleon and 
A magnetic moments is derived to next to next to leading nonanalytic order. Numerical 
simulation results for the proton and A baryon magnetic moments in quenched QCD 
reveal dramatic signatures of the quenched meson cloud, which are in accord with the 
^ i predictions of quenched chiral perturbation theory. 
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1. NUMERICAL SIMULATIONS WITH FLIC FERMIONS 



^ [ Access to leading-edge supercomputing resources coupled with advances in the for- 
mulation of computationally- inexpensive chirally-improved lattice fermion actions [ |2] 
enable the numerical calculation of hadron structure in the chiral regime. In this regime, 
the pseudo-Goldstone boson dressings of hadrons give rise to significant non-analytic cur- 
vature in the quark-mass dependence of observables. The magnetic moments of baryons 
have been identified [El 111 as providing an excellent opportunity for the direct observation 
of chiral nonanalytic behavior in lattice QCD, even in the quenched approximation. 

The numerical simulations of the electromagnetic form factors presented here are carried 
out using the Fat Link Irrelevant Clover (FLIC) fermion action [ |2] in which the 
irrelevant operators introduced to remove fermion doublers and lattice spacing artifacts 
are constructed with smoothed links. These links are created via APE smearing [E]; a 
process that averages a link with its nearest transverse neighbors in a gauge invariant 
manner. Iteration of the averaging process generates a "fat" link. 

The use of links in which short-distance fluctuations have been removed simplifies the 
determination of the coefficients of the improvement terms in both the action and its 
associated conserved vector current. Perturbative renormalizations are small for smeared 
links and the mean-field improved coefficients used here are sufficient to remove 0{a) 
errors, in the lattice spacing a, from the lattice fermion action. The key is that both the 
energy dimension-five Wilson and Clover terms [ E] are constructed with smooth links, 



2 



while the relevant operators, surviving in the continuum limit, are constructed with the 
original untouched links generated via standard Monte Carlo techniques. 

FLIC fermions provide a new form of nonperturbative 0{a) improvement [I21IZI where 
near-continuum results are obtained at finite lattice spacing. Access to the light quark 
mass regime is enabled by the improved chiral properties of the lattice fermion action. 
The magnitude of additive mass renormalizations is suppressed [[7j which otherwise can 
lead to singular behavior in the propagators as the quarks become light. 

The 0(a)-improved conserved vector current is used. Nonperturbative improvement 
is achieved via the FLIC procedure where the terms of the Noether current having their 
origin in the irrelevant operators of the fermion action are constructed with mean-field 
improved APE smeared links. The preliminary results presented here are from a sample 
of 255 20^ X 40 mean- field improved Luscher-Weisz [H] gauge field configurations having 
a lattice spacing of 0.128 fm as determined by the Sommer scale Tq = 0.50 fm. 

2. CHIRAL NONANALYTIC BEHAVIOR 

The truncation of the low-energy expansion of chiral effective field theory introduces 
errors into the predictions of chiral perturbation theory (^PT). In the process of a simple 
truncation, one sets the coefficients of higher-order terms of the expansion (both analytic 
and nonanalytic) to zero by hand. While such a procedure is often described as "sys- 
tematic" or "model independent" the truth is that the coefficients of these higher-order 
terms are generally not zero. Thus, the truncated expansion is a poor representation of 
the chiral expansion of QCD. In this case, "model independent" simply means that no 
attempt has been made to estimate the coefficients of higher order terms in the expansion. 
It is not necessarily a good feature. 

While such an approach might be forgiven if the coefficients of the higher order terms 
were indeed small, there is now mounting evidence that this is not the case. For the 
nucleon mass, the best determination of the low-energy constants from the physical 
nucleon mass and state-of-the-art lattice QCD results [ indicate the nucleon mass has 
the following chiral expansion (in appropriate powers of GeV) 



where Xi the known model-independent coefficients of the leading nonanalytic terms 
of the expansion and the quoted uncertainties are purely systematic [E]. The coefficient 
of the term from the Nn self energy alone is — 75 ± 35 GeV~^. Hence the 

"systematic" approach of setting the coefficients of all higher order terms to zero can be 
troublesome for the power-series like expansion of dimensional- regularization (DR). 

Fortunately there is a way to estimate the coefficients of the higher order terms, while 
preserving the model-independent features of chiral effective field theory to the chiral- 
counting order that one is working. Through the process of regulating loop integrals via 
a finite-range regulator (FRR) [CH], one re-sums the chiral expansion in a manner which 
preserves the model-independent features of chiral perturbation theory. For example, 
the coefficients of the leading nonanalytic terms of the expansion are preserved exactly. 
However, in the process of expanding the expressions of FRR xPT to recover the expansion 



TTiN = (0.897 ± 0.001) + (2.84 ± 0.04) ml + X3rnl + (22.0 ± 1.6) 




(1) 
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(c) = -(b) - -(a) (b) = (a) 

Figure 1. Quark-flow diagrams for meson-cloud contributions to A^"*" in full QCD. 

of DR in Eq. (P), one also encounters higher-order terms of the chiral expansion, whose 
coefficients are functions of the regulator parameter (A) governing the finite range of the 
regulator. By optimizing A in a fit of the FRR expansion to lattice QCD data, one obtains 
estimates for the higher order terms of the chiral expansion, while maintaining complete 
model-independence to the chiral order one is working. 

The introduction of a finite range regulator opens the question of the functional form of 
the regulator. Clearly one needs a functional form that preserves the low energy physics 
of the chiral expansion, while suppressing high-energy contributions where the internal 
structure of the effective fields becomes important. In practice, sharp-cutoff, monopole, 
dipole and Gaussian vertex regulators have been investigated in some detail. Provided 
one allows the regulator parameter to be constrained by lattice QCD data, a remarkable 
robustness is observed in the predictions of the chiral expansion. Systematic errors for the 
nucleon mass have been estimated at less than 1%, provided a smooth FRR (monopole, 
dipole, Gaussian) is selected [UHl- These regulators not only provide estimates for the 
coefficients of higher-order analytic terms of the expansion, but also provide estimates 
for the coefficients of non-analytic terms which must also appear in the expansion. This 
latter feature of the smooth regulators has been identified as key to the success of the 
smooth regulators over the sharp cutoff [HU]. 

We use the diagrammatic method for evaluating the quenched chiral coefficients of 
leading nonanalytic terms in heavy-baryon quenched xPT [ El E] • Results for the proton 
magnetic moment to next to next to leading nonanalytic (NNLNA) order [ HI ^] are 
generalized to the FRR approach used here. In quenched QCD, the A form factors are 
simply proportional to the charge of the baryon [El- Hence, consideration of the A"'"^ 
charge state is sufficient to determine the chiral expansion for all charge states. 

The presence of the A Ntt decay channel is particularly important for the mass 
dependence of A properties. Rapid curvature associated with nonanalytic behavior is 
shifted to larger pion masses near the A^-A mass splitting, m-jr ~ Ma — M^. As described 
below, quenched-QCD decay-channel contributions come with a sign opposite to that of 
full QCD. This artifact holds tremendous promise for revealing unmistakable signatures 
of the quenched meson cloud. 

The change in sign for the decay-channel contributions is easily understood through 
the consideration of the quark fiow diagrams of Fig. ^ illustrating the meson-cloud con- 
tributions to the A^^ resonance in full QCD. Quark fiow diagram (a) corresponds to 
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Figure 2. FLIC fermion simulation results for the masses of the proton (o) and A res- 
onance (A) in quenched QCD. Solid curves indicate the fits of FRR quenched xPT to 
the lattice simulation results, while dashed curves indicate the one-loop correction to the 
quenched approximation [US]. Stars denote the physical values. 

the hadronic process described at left. Since QCD is flavor-blind, the process illustrated 
in diagram (b) is equivalent to diagram (a) provided the masses of the u and d quarks 
are taken to be equal. On its own, diagram (b) describes the decay of the A"*""*" to a 
doubly-charged uuu "proton," which we denote Of course, such states do not exist 
in full QCD and diagram (c) makes contributions exactly equal but opposite in sign to 
diagram (b) when the intermediate state is a uuu proton. Upon quenching the theory, 
both diagrams (a) and (b) are eliminated, leaving only diagram (c). Hence the physics of 
the A Nti decay is present in the quenched approximation [[T3] but its contribution 
has the wrong sign. 

This aspect of the quenched A is the predominant feature giving rise to the flattening 
of the A mass as a function of quark mass at the lightest FLIC-fermion quark masses 
depicted in Fig. |21 Double-hairpin t]' contributions impact the A mass at much lighter 
quark masses most notably between the chiral and physical pion mass. 

As for the proton, the double-hairpin rj' dressing A^+ A^+r/' provides the LNA 
contribution to the A^"*" magnetic moment, generating a logarithmic divergence in the 
chiral limit. For electromagnetic form factors, Fig.^c) indicates that the coefficient of the 
NLNA contribution proportional to m,r will vanish in the quenched approximation because 
of the neutral charge of the meson. However, Fig. ^c) will make significant contributions 
when the electromagnetic current couples to the intermediate We estimate the 

tree-level magnetic moment of the uuu proton using standard SU(6) symmetry, fip++ = 
— = = |/iA++- The quark fiow diagrams of Fig. [T]also include contributions 
from A intermediate states. In terms of the full QCD process A"*""*" — ^ A++7r° we find the 
total quenched contribution to be (4/3) (A++ A++7r°). 

Figure El displays FLIC fermion simulation results for the magnetic moments of the pro- 
ton and A^ resonance in quenched QCD. The curves illustrate the fits of FRR quenched 
xPT to the lattice simulation results. Here, the analytic terms of the chiral expansion have 
been re-summed in a Fade designed to reproduce the Dirac moment mass dependence, 
eh/ 2m, at moderately large pion mass. 
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Figure 3. FLIC fermion simulation results for the magnetic moments of the proton (o) 
and A"*" resonance (A) in quenched QCD and the associated fits of FRR quenched xPT 
to the lattice simulation results. 

At large pion masses, the A moment is enhanced relative to the proton moment in 
accord with earlier lattice QCD results [ 221 1^ and model expectations. However as 
the chiral regime is approached the nonanalytic behavior of the quenched meson cloud is 
revealed, enhancing the proton and suppressing the A+ in accord with the expectations 
of quenched xPT. The quenched artifacts of the A provide an unmistakable signal for the 
onset of quenched chiral nonanalytic behavior. 
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